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Abstract 

We discuss constrained and semi-constrained versions of the next -to-minimal su- 
persymmetric extension of the Standard Model (NMSSM) in which a singlet Higgs 
superfield is added to the two doublet superfields that are present in the minimal ex- 
tension (MSSM). This leads to a richer Higgs and neutralino spectrum and allows for 
many interesting phenomena that are not present in the MSSM. In particular, light 
Higgs particles are still allowed by current constraints and could appear as decay prod- 
ucts of the heavier Higgs states, rendering their search rather difficult at the LHC. We 
propose benchmark scenarios which address the new phenomenological features, con- 
sistent with present constraints from colliders and with the dark matter relic density, 
and with (semi-)universal soft terms at the GUT scale. We present the corresponding 
spectra for the Higgs particles, their couplings to gauge bosons and fermions and their 
most important decay branching ratios. A brief survey of the search strategies for these 
states at the LHC is given. 
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1. Introduction 



The next-to-minimal supersymmetric extension of the standard model (NMSSM), in which 
the spectrum of the minimal supersymmetric extension (MSSM) is extended by one sin- 
glet superfield, has been discussed since the early days of supersymmetry (SUSY) model- 
building [1-3]. In the last decade, the NMSSM gained a renewed interest in view of its 
positive features as compared to the widely studied MSSM. Firstly, the NMSSM naturally 
solves in an elegant way the so-called fi problem [4] of the MSSM: to have an accept- 
able phenomenology, a value in the vicinity of the electroweak or SUSY breaking scale is 
needed for the supersymmetric Higgs mass parameter n; this is automatically achieved in 
the NMSSM, since the /i-parameter is dynamically generated when the singlet Higgs field 
acquires a vacuum expectation value of the order of the SUSY breaking scale, leading to a 
fundamentalLagrangian that contains no dimensionful parameters apart from the soft SUSY 
breaking terms. Secondly, as compared to the MSSM, the NMSSM can induce a richer 
phenomenology in the Higgs and neutralino sectors, both in collider and dark matter (DM) 
experiments: on the one hand, heavier Higgs states can decay into lighter ones with siz- 
able rates [5-16] and, on the other hand, a new possibility appears for achieving the correct 
cosmological relic density [17] through the so-called "singlino", i.e. the fifth neutralino of 
the model, which can have weaker-than-usual couplings to standard model (SM) particles. 
Thirdly, the NMSSM needs somewhat less fine tuning [8, 18] (although some fine tuning is 
still required [19]): the upper limit on the mass of the lightest CP-even Higgs particle is 
larger than in the MSSM, and therefore more SUSY parameter space survives the bounds 
imposed by the negative Higgs boson searches at LEP; furthermore, possible unconventional 
decays of the SM-like Higgs scalar allow it to have a relatively small mass, well below the 
SM Higgs mass limit of 114 GeV, still consistent with LEP constraints. 

Given the possibly quite different phenomenology in the Higgs sector as compared to 
the SM and the MSSM [20], it is important to address the question whether such NMSSM 
specific scenarios will be probed at the LHC In particular, it would be important to extend 
the so-called "no- lose theorem" of the MSSM [21], which states that at least one MSSM 
Higgs particle should be observed at the LHC for the planned integrated luminosity, to 
the NMSSM [6, 22] or try to define regions of the NMSSM parameter space where more 
Higgs states are visible at the LHC than those available within the MSSM [23]. However, 
a potential drawback of the NMSSM, at least in its non constrained versions, is that it 
leads to a larger number of input parameters in the Higgs sector to deal with, compared to 
the MSSM. In particular, it is clearly unfeasible to perform multi-dimensional "continuous" 
scans over the free inputs of the NMSSM, especially if each sampled point is subject to a 
complete simulation in order to be as close as possible to the experimental conditions. 

An alternative approach, acknowledged by both the theoretical and experimental com- 
munities, is that of resorting to the definition of so-called "benchmark points" (or slopes, 
or surfaces) in the SUSY parameter space [24]. These consist of a few "discrete" parame- 
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ter configurations of a given SUSY model, wfiicli embody tlie most peculiar/representative 
pfienomenological features of the model's parameter space. Using discrete points avoids 
scanning the entire parameter space, focusing instead on representative choices that reflect 
the new interesting features of the model, such as new signals, peculiar mass spectra, etc. 
A reduced number of points can then be subject to full experimental investigation, without 
loss of substantial theoretical information. 

While several such benchmark scenarios have been devised for the MSSM [24, 25] and 
thoroughly studied in both the collider and the DM contexts, limited progress has been made 
so far in this direction for the case of the NMSSM. In Ref. [9], an earlier attempt was made 
to address the possibility of establishing a no-lose theorem for the NMSSM Higgs sector 
at the LHC through benchmark points (in a non-universal low energy setup in which the 
SUSY particles are very heavy). However, the corresponding spectrum was not consistent 
with DM constraints (in particular on the lightest SUSY particle, LSP), which were not 
yet established for the NMSSM at the time. In addition, many of the points discussed in 
Ref. [9] have become ruled out due to the new lower value of the top quark mass, and to 
more stringent constraints from collider searches and precision measurements. Also the tools 
to calculate the Higgs and SUSY particle spectra have been upgraded since then [26,27]. 

In this report, we build on the experience of Ref. [9] and define benchmark points which 
fulfill the present collider and cosmological constraints using the most up-to-date tools to 
calculate the particle spectra. However, in contrast to Ref. [9] , we work in the framework of 
a (semi-) constrained NMSSM parameter space, henceforth called cNMSSM, where the soft 
SUSY-breaking parameters are defined at the Grand Unified (GUT) scale: on the one hand, 
this approach leads to a much more plausible sparticle spectrum and allows one to relate 
features of the Higgs sector to properties of the neutralino sector (and hence of the LSP); on 
the other hand, this restricted parameter space still contains the phenomenological features 
of the NMSSM that are very distinctive from those of the MSSM, and is suitable for inten- 
sive phenomenological investigation by the experimental collaborations. The emphasis will 
primarily be on the different possible scenarios within the Higgs sector and the implication 
for Higgs searches at the LHC. However, we will also comment on the possible implications 
of these benchmark points for the cosmological relic density of the lightest neutralino DM 
candidate. Finally, we will describe the tools used to define such benchmark scenarios. 

The report is organized as follows. In the next section, we define the NMSSM and its 
particle content with some emphasis on constrained scenarios defined in terms of soft terms 
specified at the GUT scale, discuss the tools that allow to calculate the particle spectrum, 
and the various constraints that should be imposed on the latter. In section 3, we propose 
five benchmark points which lead to Higgs sectors that are different from those of the MSSM 
and discuss their main features. In section 4, we outline the possible search strategies at the 
LHC for the Higgs particles in these scenarios. A brief outlook is given in section 5. 
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2. The NMSSM and its particle spectrum 
2.1 The unconstrained NMSSM 

In this paper, we confine ourselves to the NMSSM with a scale invariant superpotential. Al- 
ternative generalizations of the MSSM - known as the minimal non-minimal supersymmetric 
SM (MNSSM), new minimally-extended supersymmetric SM or nearly-minimal supersym- 
metric SM (nMSSM) or with additional U(l)' gauge symmetries - exist [28], but these will 
not be considered here, nor the case of explicit CP violation [29]. The scale invariant super- 
potential of the NMSSM is given, in terms of (hatted) superfields, by 

W = + I + KQHS'r - hQHab'ji - KLH^r^^ (1) 

in which only the third generation fermions have been included (with possible neutrino 
Yukawa couplings have been set to zero), and Q, L stand for superfields associated with the 
{t, h) and (r, z/^) SU(2) doublets. The first two terms substitute the iiHuHd term in the 
MSSM superpotential, while the three last terms are the usual generalization of the Yukawa 
interactions. The soft SUSY breaking terms consist of the scalar mass terms for the Higgs 
and sfermion scalar fields which, in terms of the fields corresponding to the complex scalar 
components of the superfields, are given by 

and the trilinear interactions between the sfermion and Higgs fields, 

- £tHi = XAxH^HdS + ^kA^S^ + htAQHJl - hhAhQHj>% - KArLRdf^ + h.c. . (3) 

In an unconstrained NMSSM with non-universal soft terms at the GUT scale, the three 
SUSY breaking masses squared for Hd and S appearing in £mass can be expressed in 
terms of their vevs through the three minimization conditions of the scalar potential. Thus, 
in contrast to the MSSM (where one has only two free parameters at the tree level, generally 
chosen to be the ratio of Higgs vacuum expectation values (vevs) tan (3 and the mass of the 
pseudoscalar Higgs boson), the Higgs sector of the NMSSM is described by the six parameters 

A , ft , , A^, tan/3 = {H,,)/{Hd) and ^efr = A(5) . (4) 

One can choose sign conventions such that the parameters A and tan/3 are positive, while 
the parameters /t. Ax-, A,^ and fi^s can have both signs. 

In addition to the above parameters of the Higgs sector, one needs to specify the soft 
SUSY breaking mass terms in eq. (2) for the scalars, the trilinear couplings in eq. (3) as well 
as the gaugino soft SUSY breaking mass parameters to describe the model completely, 
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Mi55 + M2^iyWa + M3^G"Ga + h.c. . (5) 
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Clearly, in the limit A ^ with finite /ieg, the NMSSM turns into the MSSM with a 
decoupled singlet sector. Whereas the phenomenology of the NMSSM for A ^ could still 
differ somewhat from the MSSM in the case where the lightest SUSY particle is the singlino 
(and hence with the possibility of a long lived next-to-lightest SUSY particle [30]), we will 
not consider this situation here. 

2.2 The constrained NMSSM 

As the number of input parameters is rather large, one can attempt to define a constrained 
model, hereafter called the cNMSSM, in which the soft SUSY-breaking parameters are fixed 
at the GUT scale, leading to only a handful of free input parameters to cope with. This 
approach is motivated by the fact that in many schemes for SUSY-breaking, the soft SUSY 
breaking parameters are predicted to be universal at a very high energy scale. For example, 
in the cMSSM or mSUGRA scenario, one imposes a common gaugino mass M1/2, a scalar 
mass mo and a trilinear coupling Aq at Mqut, leading to only four continuous free parameters 
(together with tan /3) and the sign of /i. The values of the numerous soft SUSY-breaking 
parameters at low energies are then obtained through the renormalization group evolution 
(RGE). Analogously, the cNMSSM allows to describe the entire sparticle spectrum, including 
the chargino and neutralino sectors, in terms of a small number of free parameters. This is 
in contrast to the usual procedure where one postulates universal sfermion masses directly 
at the weak scale which seems less plausible since these masses would then be non-universal 
at the high (GUT) scale. 

In analogy to the cMSSM, one can impose unification of the soft SUSY-breaking gaugino 
masses, sfermion and Higgs masses and trilinear couplings at the scale Mqut^ 

Mi,2,3 = Ml/2 , rnp^= rnu, = mo , Aj = Aq . (6) 

The cMSSM has two additional parameters (/i and B) beyond those in eq. ([6]), but the 
correct value of M^o imposes one constraint (typically used to determine |/i|) and B can be 
replaced by tan/5. Likewise, the fully constrained NMSSM has two additional parameters 
(A and k) beyond the three parameters in eq. ([6]), and the correct value of M^o imposes 
one constraint. Hence, a priori the number of free parameters in the cMSSM and the fully 
constrained NMSSM is the same. 

In principle, one could minimize the effective potential of the cNMSSM with respect to 
Hu, Hd and S", and determine the overall scale of the soft terms in eq. ([6]) from the correct 
value of M^o: this approach has been pursued in Ref. [2]. However, since tan/3 is then 
obtained as output (while the top quark Yukawa coupling ht is an input), it becomes very 
difficult to obtain the correct value for mtop- Also, the numerical minimization of the effective 
potential including radiative corrections is quite computer-time consuming. 

Therefore it is much more convenient to employ the analytic form of the three mini- 
mization conditions of the effective potential of the NMSSM: for given M^o, tan/3, A and 
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all soft terms at the weak scale (except for m|), these can be solved for \fica\ (or | {S) |), 
K and m|. As in the cMSSM, the sign of fi^s can still be chosen at will; in some sense, 
the determination of the parameter B through the minimization conditions of the MSSM is 
replaced by the determination of k in the NMSSM and the additional minimization w.r.t. 
(S) in the NMSSM leads to the determination of m|. 

Assuming the weak scale soft scalar masses (apart from m|) to arise from a unified scalar 
mass TTiQ, the determined weak scale value of m| - when run up to the GUT scale - will not 
coincide with ttiq in general. One could confine oneself to regions in parameter space where 
the difference between m| and rriQ is negligibly small, but we found that the phenomenology 
of this fully constrained NMSSM hardly differs (at least in the Higgs sector on which we 
are focusing here) from the one of the cMSSM once present LEP constraints are imposed. 
Hence we will relax the hypothesis of complete unification of the soft terms in the singlet 
sector since the singlet could play a special rol^, and we will allow for both m| 7^ itiq and 

7^ Aq at Mgut- Note that, although m| at Mqut can be negative, this does not signal 
an instability of the potential: the direction (S) — 00 is always protected by a quartic self 
coupling ~ k"^, which leads to (S) ~ Msusy at the minimum of the potential. 

In addition, for some of the benchmark points (see points P4 and P5 below), we will also 
relax the unification hypothesis for mjj^ and rn'jj^, in analogy to corresponding scenarios 
within the MSSM (called non-universal Higgs model or NUHM [32]). The hypothesis Ax = 
Aq will also be relaxed for point 5. Such points in parameter space can have additional 
unconventional properties, whose phenomenology should also be investigated. 

2.3 The Higgs and SUSY spectra 

Following the procedure employed by the routine NMSPEC within NMSSMTools [27], which 
calculates the spectra of the Higgs and SUSY particles in the NMSSM, a point in the 
parameter space of the cNMSSM is defined by the soft SUSY breaking terms at Mgut (except 
for the parameter m|) , tan f] at the weak scale, A at the SUSY scale (defined as an average 
of the first generation squark masses) and the sign of the parameter fics- The parameters k, 
m| and \nes\ are determined at the SUSY scale in terms of the other parameters through 
the minimization equations of the scalar potential. 

The RGEs for the gauge and Yukawa couplings have to be integrated from the weak 
scale up to Mqut (defined by the unification of the gauge couplings gi and (72)5 and the 
RGEs for the soft terms from Mqut down to the weak scale. Since k and m| are computed 
at the weak scale, and threshold effects for the gauge and Yukawa couplings depend on 
the soft terms at the SUSY scale, some iterations are necessary in order to satisfy the 
desired boundary conditions for the soft terms at the GUT scale, but usually the procedure 
converges quite rapidly (at least for not too large values of tan/5 or A). For the most 

^An example is when the singlet is mixed with or identified with a radion originating from a 5d brane 
world [31]; however, here we do not consider a Kaluza-Klein mass scale below A/qut- 
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relevant SM parameters, the strong coupling and the top/bottom quark masses, we chose 
a,{Mzo) = 0.1172, mb{mb)^ = 4.214 GeV and m^°^ = 171.4 GeV. 

After RGE running is completed, the gluino, chargino, neutralino and sfermion masses 
are computed including dominant one loop corrections to their pole masses. The lightest 
scalar Higgs pole mass is determined to the following accuracy (in a notation where its 
tree level mass is proportional to 0{g'^), where g denotes any of the electroweak gauge 
couplings): one loop corrections of 0{hffj) and Oihl^^g"^) are computed exactly; one loop 
corrections of the orders g^, fi'^A^, A^, (^^k^, and X^n^ include only terms involving large 
logarithms ln(Mj^/M|o), where Mi are potentially large Higgs or sparticle masses, while two 
loop corrections of 0{h^f,) and 0{hjf,as) include only terms with two powers of large logs. 

Once the pole masses are known, all Higgs decay branching ratios into SM and SUSY 
particles (an adaptation of the decays in the MSSM [33] to the NMSSM case) are determined 
including dominant (mainly QCD) radiative corrections [34] , as well as sparticle loops which 
contribute to the couplings of a neutral Higgs to two photons or gluons [35]. 

When the spectrum and the couplings of the Higgs and SUSY particles are computed, 
available Tevatron and LEP constraints are checked. The results of the four LEP collabora- 
tions, combined by the LEP Higgs working group, are included [36]. More specifically, the 
following experimental constraints are taken into account: 

(i) The masses of the neutralino as well as their couplings to the boson are compared 
with the LEP constraints from direct searches and from the invisible boson width; 

(ii) Direct bounds from LEP and Tevatron on the masses of the charged particles {h^, x^, 
q, I) and the gluino are taken into account; 

(iii) Constraints on the Higgs production rates from all channels studied at LEP. These 
include in particular Z^h^ production, being any of the CP-even Higgs particles, 
with all possible two body decay modes of (into b quarks, r leptons, jets, photons 
or invisible), and all possible decay modes of /i° of the form '^i^i' ^'j being 
any of the CP-odd Higgs particles, with all possible combinations of decays into 
b quarks, c quarks, r leptons and jets. Also considered is the associated production 
mode e"'"e~ — > /i°a° with, possibly, — ^ '^j'^j'- (I^ practice, for our purposes, only 
combinations of z = 1, 2 with j = 1 are phenomenologically relevant.) 

We stress that light Higgs scalars (with M^^o^ ^114 GeV) can still be allowed by LEP 
constraints, if either (i) the Z^-Z^-h\ 2 coupling is heavily suppressed (if, for instance, the 
state is dominantly a gauge singlet); or (ii) M„o < 11 GeV such that the state /i5,2 decays 
dominantly into states, but the bb decay of a? is impossible; constraints from the decays 
^ CLidi — > 4r allow for M/jO ^ down to ~ 86 GeV. It is important to note, however, that 
LEP constraints are implemented only for individual processes, and that combinations of 
different processes could potentially rule out seemingly viable scenarios. 

Finally, experimental constraints from B physics [37] such as the branching ratios of the 
rare decays BR(i? Xs^), BR(i?s /U"*"/^") and BR(5"'" t^i^t) and the mass differences 
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AMs and AM^, are also implemented; compatibility of each point in parameter space with 
the current experimental bounds is required at the two sigma level. 

The new features of the NMSSM have an impact on the properties of the lightest neu- 
tralino as a dark matter candidate [17]. For instance, given the presence of a fifth neutralino 
(singlino), the composition of the annihilating WIMPs can be significantly different from 
those within the MSSM in wide regions of the parameter space. In particular, the LSP 
can have a large singlino component, in which case one has new couplings of the LSP to 
singlet-like Higgs states whose mass can be substantially lighter than the Higgs states within 
the MSSM. In the presence of light and states, there are new channels through which 
neutralino annihilation can occur in the NMSSM: hi, h\h\, h\a\ and a^a^, either via 
s-channel Z^, a°, or t-channel neutralino exchange [17]. Among the proposed benchmark 
points, we will find examples of these new features. 

Technically, the relic abundance of the NMSSM dark matter candidate x\ is evaluated 
via a link to MicrOMEGAS [38]. All the relevant cross sections for the lightest neutralino 
annihilation and co-annihilation are computed. The density evolution equation is numerically 
solved and the relic density of Xi is calculated. The result is compared with the "WMAP" 
constraint 0.094 < VtcBuh'^ < 0.136 at the 2a level [39]. 

3. The benchmark points 

As already mentioned in the introduction, in view of the upcoming LHC, quite some work 
has been dedicated to probe the NMSSM Higgs sector over the recent years. In the NMSSM, 
two different types of scenarios have been pointed out, depending on whether Higgs-to-Higgs 
decays are kinematically allowed or forbidden. Within the first category, where Higgs-to- 
Higgs decays are kinematically allowed and can be dominant, there are two possibilities, 
each associated with light scalar /pseudoscalar Higgs states: 

(i) The hghtest pseudoscalar Higgs boson a? is rather light, M^o < 40-50 GeV, and the 
lightest CP-even Higgs particle h\ has enough phase space for the decay into two pseudo- 
scalar Higgs particles, h\ — ^ a^a?, to be kinematically accessible. In this case, the branching 
ratio for the decay h\ a^al is typically very large and this new decay channel is the dom- 
inant one. Concerning the lightest pseudoscalar decays, there are two further possibilities: 
either M^o > 10 GeV and the boson decays into a pair of b quarks or a pair of tau leptons 
(with the former decay being in general dominant), leading to hi — >• a^a? 4r, 46 and 2r2fe 
final states, or M^o < 10 GeV and the dominant decay mode of the boson is into a pair 
of tau leptons, leading to hi a^a? ~^ 4r final states. In this latter scenario, one can still 
distinguish two situations for the h^ boson: its mass is either close to its theoretical upper 
limit of 130 GeV (within the cNMSSM), or to the lower limit of 90 GeV, still in agreement 
with constraints from Higgs boson searches at LEP2. 

(ii) The lightest CP-even Higgs boson is relatively light, Mj^o < 50 GeV, and decays into 
bb pairs. (The situation where Mj^o < 10 GeV, such that the latter channel is kinematically 



8 



closed and the decay t^t~ dominates, is very constrained by LEP data.) In this case, 
the next-to-hghtest CP even Higgs boson h\ is SM-hke, and has a mass below ~ 150 GeV. 
Even so, it can still decay into two h\ bosons leading to the final topologies — > h\h\ 4r, 
2r26 and 46, the latter final state being largely dominating. 

The second category of scenarios, where Higgs-to-Higgs decays are suppressed, includes 
regions of the parameter space where all Higgs particles are relatively light with masses in 
the range 90-200 GeV. This opens the possibility of producing all the five neutral and the 
charged Higgs bosons at the LHC This scenario is similar to the so-called "intense coupling 
regime" of the MSSM [40,41] (which has been shown to be rather difficult to be fully covered 
at the LHC), but with two more neutral Higgs particles. 

In the context of the general NMSSM (without unification constraints), there exists 
also a decoupling regime in which all the Higgs bosons are heavy and decouple from the 
spectrum except for one CP-even Higgs particle, whose mass can be up to ~ 140 GeV [42]. 
Such decoupling scenarios are difficult to realize within the cNMSSM and, in any case, they 
would be very difficult to disentangle from an MSSM "M™*™" -scenario at the LHC, since 
the separation of the two would require a determination of the parameter tan/5 which is 
notoriously difficult to achieve at the LHC. 

In the following, we propose five benchmark points of the NMSSM parameter space, PI 
to P5, in which the above discussed scenarios are realized. Each point is representative 
of distinctive NMSSM features: points PI to P3 exemplify scenarios where h\ decays into 
light pseudoscalar states (decaying, in turn, into hh or t^t~ final states), P4 illustrates the 
NMSSM possibility of a very light while P5 corresponds to the case where all Higgs bosons 
are rather light. In all cases, the input parameters as well as the resulting Higgs masses and 
the relevant decay information are given in Tabled! In the subsequent Table [2], we summarize 
the most relevant features of each point regarding the cosmological relic density. 

As discussed in Ref. [27], the scenarios where the pseudoscalar Higgs boson is light and 
the decay h\ a^a^ is dominant, can be realized within the cNMSSM with nearly universal 
soft terms at the GUT scale, the exception being the parameters m| and which are 
chosen independently. A corresponding region in parameter space, containing PI, P2 and 
P3, is shown in Fig. [H where we chose tuq = 174 GeV, M1/2 = 500 GeV, Aq = -1500 GeV, 
tan/3 = 10 and the sign of fics is positive. As a function of the parameter A, we show the 
values of Af^ that are allowed by LEP constraints, the masses of the and a? bosons (where 
larger values of M^o correspond to larger values of the trilinear coupling \Af^\), and the dark 
matter relic density QcDuh"^- The benchmark points PI, P2 and P3 discussed below are 
chosen at the upper and lower boundaries of this region of the cNMSSM parameter space. 
The lower limit on the pseudoscalar mass, M^o > 8 GeV, in the last plot of Fig. 1 originates 
from constraints from B physics; see Ref. [37] for details. 

In Fig. [2], we show the obtained results from a scan in the [A, A^] parameter space (with 
the values of the other parameters fixed to those used in Fig. 1) for the masses of the CP- 



9 




-300 
-400 



.1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I- 

;P1,P2 I I I I : 



.1 0.15 0.2 0.25 0.3 0.35 0.4 
X 

.P3 




I I I I I I I I I I I I I I I I I I I I I 



0.1 0.15 0.2 0.25 0.3 0.35 



D. h 



0.095^ 



I I I I I I I I 1^ 




P1,P2 

•I I I I I I I I I I I I I I I I I I I I I I I I I I I I I- 



M 



a[' 30 

20 



1 0.15 0.2 0.25 0.3 0.35 0.4 



P2b-'. 

Q Fi I I I !■ I- 1 I > I I I I I I I I I I I I I I I I I I I 



0.1 0.15 0.2 0.25 0.3 0.35 
X 



0.4 



Figure 1: Allowed values of A^, Mf^o and M^o (in GeV) as well as f2cDM^^ as a function 
of A. We take mo = 174 GeV, M1/2 = 500 GeV, Aq = -1500 GeV and tan/3 = 10; m| is 
determined from the electroweak symmetry breaking conditions. 

even and CP-odd states, as well as of the branching ratios of the decays /i^ a^a^ 
and tti T^T^ . In the two upper frames for the Higgs masses, a scan is performed over the 
same [A, A^] range as in Fig. [H which therefore includes the three points PI, P2 and P3. In 
the lower frames for the Higgs masses and for the branching ratios, a finer scan zooms on 
the [A, A^] range which involves only the two points PI and P2. 

A very light state (and, thus, with the decay /12 — ^ /i^/i? being the dominant channel) 
is represented by P4, which can be obtained once one relaxes the universality conditions on 
the soft SUSY-breaking Higgs mass terms, 7^ i^Hu 7^ '^o- A scenario as in P5, in which 
all NMSSM Higgs bosons are light, is possible if one allows additionally for Ax 7^ Aq. 

For all benchmark points the numerical value of the pseudoscalar Higgs boson mass M^o is 
quite sensitive not only to the NMSSM input parameters (notably to the trilinear coupling 
Ak), but also to the employed SM parameters and to the precision with which radiative 
corrections are computed. Since the numerical values of M^o are phenomenologically much 
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Figure 2: Scans in the [A, A^] parameter space for the h^, masses, and the branching ratios 
for the decay modes a^a^ and a? t^t^ . The two upper frames include PI to P3 

while the remaining frames include only PI and P2; other parameters are as in Fig. 1. 

more important than, for instance, A,^ at Mqut, we consider the benchmark points to be 
defined in terms of M^o rather than in terms of at Mqut- Next, we summarize the most 
relevant phenomenological properties of the benchmark points. 
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In the first two points PI and P2, the lightest CP-even Higgs boson has a mass of 
Mi^o ^120 GeV and is SM-hke, as reflected by the corresponding couplings to gauge bosons 
Ri, top quarks ti and bottom quarks bi, which are almost equal to unity, when normalized 
to the SM Higgs boson couplings: see Table [TJ The lightest CP-odd Higgs boson has masses 
of 40.5 GeV and 9.09 GeV, which are obtained by choosing at Mqut as —33.9 GeV 
and —33.4 GeV, respectively; the remaining parameters of the Higgs sector at Msusy, like 
K, Ax, fics and the gaugino mass parameter M2, are given in the second part of Tab. [TJ 

In both cases PI and P2, the decay channel a^a? is largely dominating with a 

branching rate very close to 90%, while the decays bb and r"*"r~ are suppressed by 

an order of magnitude as compared to the SM. The most relevant difference between the 
two scenarios concerns the mass and decays of the lightest pseudoscalar state. In PI the 
pseudoscalar a° decays into b quarks and r leptons, with branching fractions of ~ 90% and 
~ 10%, respectively. On the other hand, in P2 the pseudoscalar with its mass M^o ~9.09 
GeV decays dominantly into r+r" pairs, with a rate close to 90%. 

For point P3, the same inputs of points PI and P2 are chosen except for the and 
A parameters, which are now varied as to have a lighter state. This again leads to 
a pseudoscalar Higgs boson which has approximately the same mass as in scenario P2, 
MaO~9.96 GeV, and which decays almost exclusively into r+r" final states. The difference 
between P3 and P2 is the lightest CP-even Higgs boson h^, which has a mass M^^o ~ 90 GeV, 
lower than in scenarios PI and P2. In this case, and although is still SM-like, exhibiting 
couplings to gauge bosons, top and bottom quarks that are very close to those of the SM 
Higgs boson, it decays nevertheless almost exclusively into pairs, with a rate close to 
100%. Another difference between P2 and P3 is that in the former case the decay mode 
hi a^Z^ is kinematically possible but the branching ratio for this new interesting channel 
that is not listed in the table is rather small, BR(/i° — > a^Z^) ~ 3%. 

Note that in all these first three points, the heaviest neutral Higgs particles /ig) ^3 and Og, 
as well as the charged Higgs states h^, all have masses close to or above 1 TeV. The main 
decay modes are into bb and ti for the neutral and tb for the charged states, as tan P is not 
too large and the tt-Higgs couplings are not very strongly suppressed, while the branching 
fractions for the neutral Higgs-to-Higgs decays, in particular the channels — >■ h^hi and 
^2 ~^ o^a?, are very tiny, not exceeding the permille level. 

Regarding the properties of the DM candidate, PI, P2 and P3 exhibit a lightest neutralino 
which is bino-like, and whose mass m^o ^208 GeV. In all three cases, the correct cosmological 
density, QcDuh"^ — 0.1, is achieved through the co-annihilation with the fi slepton, which 
has a mass comparable to that of the LSP, m^^ ^ 213, 213 and 215 GeV, respectively, for 
PI, P2 and P3. In each benchmark scenario, the dominant co-annihilation channels are thus 
X? h IT (~ 33%) and x? h ^ (~ 10%). 

The benchmark point P4 corresponds to a scenario in which the CP-even boson h\ is 
very light, M^jO = 32.3 GeV, and singlet-like, with a singlet component above 99%. The 
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Table 1: Input and output parameters for the five benchmark NMSSM points. 



Point 


PI 


P2 


P3 


P4 


P5 






sisnf iinfr) 

'-'^o^^V( err J 


+ 


+ 


+ 




+ 


tan /3 


10 


10 


10 


2.6 


6 


mo (GeV) 


174 


174 


174 


775 


1500 


A/i/2 (GcV) 


500 


500 


500 


760 


175 


An 


-1500 


-1500 


-1500 


-2300 


-2468 


Ax 


-1500 


-1500 


-1500 


-2300 


-800 


A^ 


-33.9 


-33.4 


-628.56 


-1170 


60 


NUHM: Mh^ (GcV) 








880 


-311 


NUHM: Mh^ (GcV) 








2195 


1910 


Parameters at the SUSY scale 




A (input parameter) 


0.1 


0.1 


0.4 


0.53 


0.016 




0.11 


0.11 


0.31 


0.12 


-0.0029 


Ax (GcV) 


-982 


-982 


-629 


-510 


45.8 


A^ (GcV) 


-1.63 


-1.14 


-11.4 


220 


60.2 


A/2 (GcV) 


392 


392 


393 


603 


140 


Mcff (GcV) 


968 


968 


936 


-193 


303 


CP even Higgs bosons 




m/jO (GcV) 


120.2 


120.2 


89.9 


32.3 


90.7 


i?i 


1.00 


1.00 


0.998 


0.034 


-0.314 


ti 


1.00 


1.00 


0.999 


0.082 


-0.305 


61 


1.018 


1.018 


0.975 


-0.291 


-0.644 


BR(/iy ^ fefe) 


0.072 


0.056 


7 X 10-* 


0.918 


0.895 


BR(/l^ ^ T+T-) 


0.008 


0.006 


7 X 10-^ 


0.073 


0.088 


BR(/i^ ^ a>^) 


0.897 


0.921 


0.999 


0.0 


0.0 


Wfto (GeV) 


998 


998 


964 


123 


118 


i?2 


-0.0018 


-0.0018 


0.005 


0.999 


0.927 


t2 


-0.102 


-0.102 


-0.095 


0.994 


0.894 


62 


10.00 


10.00 


9.99 


1.038 


2.111 


BR(/i^ ^ hh) 


0.31 


0.31 


0.14 


0.081 


0.87 


BR(/7,^ tt) 


0.11 


0.11 


0.046 


0.0 


0.0 


rSr(,^Ai2 — > fljZ j 


0.23 


0.23 


0.72 


0.0 


0.0 




2142 


2142 


1434 


547 


174 


CP odd Higgs bosons 




TO^o (GcV) 


40.5 


9.1 


9.1 


185 


99.6 




0.0053 


0.0053 


0.0142 


0.0513 


-0.00438 




0.529 


0.528 


1.425 


0.347 


-0.158 


BR(a^ bb) 


0.91 


0. 


0. 


0.62 


0.91 


BRK ^ T+r-) 


0.085 


0.88 


0.88 


0.070 


0.090 


niao (GcV) 


1003 


1003 


996 


546 


170 


Charged Higgs boson 




mh± (GcV) 


1005 


1005 


987 


541 


188 
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Table 2: LSP properties and relic density for the five benchmark NMSSM points. 



Point 


PI 


P2 


P3 


P4 


P5 


Dark matter 




LSP (xi) mass (GcV) 


208 


208 


208 


101 


70.4 


Nn 


0.999 


0.999 


0.999 


-0.039 


0.977 


Nl2 


-0.008 


-0.008 


-0.009 


0.043 


-0.098 




0.048 


0.048 


0.050 


-0.028 


0.178 


Nii 


-0.015 


-0.015 


-0.016 


0.405 


-0.068 


Ni5 








0.003 


0.912 


-0.003 




0.099 


0.099 


0.130 


0.099 


0.105 



lightest scalar Higgs particle predominantly decays into bb pairs, with BR(/i5 bb) = 92%, 
and to a smaller extent into r pairs with BR(/i'j' — > r+r^) ^ 7%. The CP-even /ig boson 
has a mass of M^^o ~ 123 GeV and is SM-like, with normalized couplings to W jZ^ and tjb 
states close to unity. However, since M^^o > 2Mi^o, it mostly decays into two bosons, with 
BR{hl h\h\) ~ 88%, and the dominant SM-like bb decay mode occurs only at a rate 
smaller than 10%. The lightest CP-odd particle is not very heavy, M^o = 185 GeV, and 
decays mostly into fermion pairs BR(a'j' — > hb) ~ 61% and BR(a5 T^r~) ~ 7%. The 
other dominant decay is the interesting channel h\Z^ which has a branching ratio of 

the order of 30%. Finally, the heaviest CP even and CP-odd a\ states and the charged 
h"^ particles have masses in the 500 GeV range and will mostly decay, as tan (3 is small, into 
tt/tb final states for the neutral/charged states. All these features make the phenomenology 
of point P4 rather different from that of points PI to P3 discussed above. 

To achieve a correct cosmological relic density, the common sfermion and gaugino mass 
parameter mo and M1/2 at the GUT scale are set close to 1 TeV. At the SUSY scale, one 
thus finds a higgsino-singlino-like neutralino LSP, whose mass is m^o ~ 100 GeV. LSP 
annihilation essentially occurs via two channels: XiXi ~^ W^W~ (60%) and XiXi ~^ Z^h\ 
(20%), mediated by s-channel Z^ and Higgs boson exchange. 

Finally, point P5 is characterized by having all Higgs particles relatively light with masses 
in the range 90 to 190 GeV. Here, the small value for the coupling A leads to a small value k ~ 
—0.003. The three CP-even Higgs bosons with masses of 91, 118 and 174 GeV, respectively, 
share the couplings of the SM Higgs boson to SM gauge bosons with the dominant component 
being taken by the /ig state. The reduced couplings of the state for this point, not given 
in the Table 1, are = —0.205, ^3 = —0.37 and 63 = 5.7. The pseudoscalar Higgs bosons 
have masses M^o ~ 100 GeV and M„o ~ 170 GeV, while the charged Higgs particle is the 
heaviest one with a mass Mh± — 188 GeV. Here, all the neutral Higgs-to-Higgs decays are 
kinematically forbidden. This is also the case of neutral Higgs decays into lighter Higgs states 
with opposite parity and gauge bosons. The only non-fermionic two-body Higgs decays are 
thus /i^ — s> Wh^ and /13 — > WW, but as the involved Higgs-gauge boson couplings are small, 
the branching ratios are tiny. 
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In this last point P5, the LSP with a mass m^o ~ 70 GeV, is a bino-like neutralino but has 
a small non-negligible higgsino component. The correct cosmological relic density ^cdm^^ — 
0.1 is achieved through the exchange of Higgs bosons in the s-channel in the annihilation 
processes XiXi ~^ bb,T~^T~ which occur at a rate of ~ 90% and 10%, respectively, i.e. 
comparable to the Higgs decay branching ratios. 

More detailed properties of the benchmark points such as more branching ratios of the 
Higgs scalars are available on the web site of Ref. [27]. 

4. Implications for the LHC 

4.1 Dominant production processes 

In the cases discussed here, at least one CP-even Higgs particle h'^ has strong enough cou- 
plings to massive gauge bosons or top/bottom quarks, i.e. Ri,ti/bi ~ 1, to allow for the 
production at the LHC in one of the main channels which are advocated for the search of 
the SM Higgs particle [20]: 

i) gluon-gluon fusion, gg — > h^, occurring through a loop of heavy top quarks which 
couple strongly to the Higgs boson, 

ii) vector boson fusion, qq qqW*W*, qqZ^*Z^* — > qqh^, which leads to two forward jets 
and a centrally decaying Higgs boson, 

in) the Higgs-strahlung processes, qq' — >• Wh^ and qq — > which lead to a Higgs and 

a massive gauge boson in the final state, 

iv) associated production with heavy quark pairs qq/ gg — * QQh'^, with Q = t,b. 

In the scenarios PI to P3, this CP-even particle is the boson which has i?i ~ ti ~ 
6i ~ 1, but which decays predominantly into a pair of light pseudoscalar Higgs particles, 
hi aja?, which subsequently decay into light fermion pairs, a? — >■ bb and r+r". In 
scenario P4, this particle is the boson which decays most of the time into a pair of 
particles, /ig — ^ ^1^?, which again decay into light fermion pairs. In these four cases, the 
backgrounds for the modes gg ^ ^ 4/ and qq/gg tth^ — > tt + Af with f = b,T, 
will be extremely large and only the vector boson fusion (owing to the forward jet tagging) 
and possibly the Higgs-strahlung (due to the leptons coming from the decays of the gauge 
bosons) can probably be viable at the LHC. Note that in the case of scenario P4, one can 
try taking advantage of the relatively light a? and the large rate for the interesting and clean 
decay mode — > h^Z^, which occurs at the level of 2% if the Z^ boson decays into electrons 
and muons. However, the cross sections for a? production in the gg — > a? ^^"^ 99 /qQ ~^ tt^^i 
processes, which are the dominant ones for the production of the MSSM CP-odd Higgs 
boson [20], are very small as the reduced a^tt coupling is very tiny, t[ ~ 0.05. One would 
have then to rely on associated production of the state with a CP-even Higgs boson, such 
as qq — ^ — > a^/i?, but the cross section is not expected to be very large. 
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In scenario P5, the particle which has couphngs to gauge bosons and top quarks that are 
close to those of the SM Higgs boson is the h2 boson which decays into bb and t~^t~ final 
states with branching ratios close to 90% and 10%, respectively. Here again, the gg fusion 
and presumably associated production with top quarks cannot be used since the interesting 
decays such as h2 —>■ WW*, Z'^Z'^* and 77 are suppressed when compared to the SM case (by 
at least a factor of four as 62 ~ 2). Thus, in this case, only the channels qq — > gg/ig qqT^T~ 
and possibly qq' Wh2 —>■ ii'bb seem feasible. The state has still non-negligible couplings 
to gauge bosons and top quarks which lead to cross sections that are "only" one order of 
magnitude smaller than in the SM. Since here, only the decays bb (90%) and t^t~ 

(10%) are again relevant, the only processes which can be used are the vector boson fusion 
and Higgs-strahlung processes discussed above. Even so, one needs a luminosity 10 times 
larger to have the same event samples as in the SM. Note that for this point, one can also 
consider associated CP-even and CP-odd Higgs production, qq Z^* h^aP^ but the cross 
sections are still small. Finally, one should note that for this point, the value of tan/3 is 
moderate and thus, the charged Higgs coupling to tb states is the smallest possible, and does 
not guarantee the detection of the bosons. 

4.2 Summary of available studies 

Let us now briefiy summarize the few detailed studies (possibly including Monte Carlo 
simulations) of the LHC potential for the NMSSM Higgs sector that have been performed 
for scenarios close to the ones discussed her^. 

In scenarios with a very light pseudoscalar Higgs boson, preliminary LHC studies focused 
on the qq — >• qqWW, qqZ^Z^ — > qqh^ qqa\a\ detection mode, i.e. via vector boson fusion 
(VBF) with forward/backward jet tagging [6]. The hope was that such NMSSM specific 
scenarios would be visible, particularly if the lightest CP-odd Higgs boson mass allowed for 
abundant a\a\ bbr^r^ decays, with both r-leptons being detected via their e,/i leptonic 
decays. At high luminosity, the VBF signal may be detectable at the LHC as a bump in the 
tail of a rapidly falling mass distribution of the 2r2j system (without b tagging). However, 
this procedure relies on the background shape to be accurately predictable. These analyses 
were based on Monte Carlo (MC) event generation via the SUSY routines of the HERWIG 
code [43] and the toy detector simulation GETJET. Further analyses based on the PYTHIA 
generator [44] and the more adequate ATLAS detector simulation ATLFAST [45] found that 
the original selection procedures may need improvement in order to extract a signal. 

Also considered was the Higgs-strahlung process with the qq' W* W^h^ — > ly^a'j'a^ 
signature, exploiting leptonic decays of gauge bosons with a subleading component from 
qq —>■ Z^* —>■ Z^h\ Z^a\a\ and, more marginally, associated production with top quarks, 
1^)99 ~^ tt^i ~^ tittia^ [12]. Regardless of the decay modes of the pseudoscalar Higgs 
bosons, it has been shown in this parton-level analysis, but in which the efficiency to trigger 

^ The prospects for the Tevatron have been discussed in Refs. [11, 13, 15] and seem quite dim. 
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on the signal is included, that at least the Higgs-strahlung process qq Wh\ should be 
taken into account along with the VBF process to improve the overall signal efficiency. This 
is particularly true for h\ masses below ~ 90 GeV where the Higgs-strahlung cross section 
exceeds that from VBF production. In another parton-level study [16], Higgs-strahlung 
with h\ — > a\a\ 46 was considered, and it was claimed that with a charged lepton for the 
W boson and the fully tagged 46 final state with Higgs mass reconstruction, the signal could 
be disentangled from the background. Nevertheless, these results need to be conffimed by 
MC and detector simulations. 

Note that in Ref. [13], the h\ —* 47 final state topology was advocated to be useful for 
the LHC at very high luminosity if the branching ratio of the decay exceeds the level of 10~^. 
However, in general, this decay has a smaller branching ratio and again, a detailed simulation 
which takes into account the experimental environment is lacking here. The scope of other 
decays, such as a^a? — >■ jjjj and jj t^t~ where j represents a light quark jet, is expected to 
be very much reduced, while the possibilities from a^a? — > 4r are currently being explored, 
in both vector boson and Higgs-strahlung production processes. 

In scenarios with a very light scalar Higgs boson such as our point P4, constraints from 
Higgs searches at LEP do not allow for h\ masses below about 10 GeV, hence the main 
decay mode would be h\ bb while the decay r^r^ would have a branching ratio of 

the order of 7-8%. The studies discussed above for very light pseudoscalar Higgs bosons but 
with the SM-like Higgs boson being the state, which then decays into two hi bosons, can 
therefore be adapted to this case. In particular, the situation for point P4 would be very 
similar to that of point PI as the production cross sections, the decay branching ratios and 
the masses of the involved primary and secondary Higgs bosons are very similar. 

Ref. [46] considered a particularly challenging NMSSM scenario, with a Higgs spectrum 
very similar to that of the MSSM, i.e., nearly degenerate (doublet dominated), heavy charged, 
scalar and pseudoscalar states and a light scalar Higgs boson at around 120-140 GeV, but 
including an additional singlet-dominated scalar and a pseudoscalar; such a scenario is some- 
what similar to our point P4. Despite of having a reasonably large production cross-sections 
at the LHC, this light Higgs boson would be difficult to observe since its main hadronic 
decays cannot be easily disentangled from the QCD backgrounds. 

In addition, some studies performed in the CP-violating MSSM, in which the boson 
can also be very light with reduced couplings to gauge bosons, can also be adapted to this 
context. Note however, that in the case of the CP-violating MSSM, one would also have a 
light boson, while in the NMSSM this state can be very heavy. 

Finally, in scenarios in which all NMSSM Higgs bosons are relatively light as in point P5, 
no detailed analysis has been performed. To our knowledge, the only study that is available 
is the one performed in Ref. [9] in which the ATLAS and CMS signal significances for the 
MSSM, assuming a high luminosity of 300 fb~^, was rescaled to take into account the re- 
duced couplings of the various Higgs particles to gauge bosons and top quarks. The effect 
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of almost overlapping resonances, new decay modes or production channels not present in 
the MSSM analyses of the ATLAS and CMS collaborations have not been considered. In 
fact, the situation in this scenario looks similar to that of the intense coupling regime of the 
MSSM [40], in which the three neutral Higgs bosons have masses in the range 100-140 GeV. 
It was shown in the detailed simulation of Ref. [41] that while it would be very difficult to 
resolve all Higgs resonances, it would be at least possible to detect one or two Higgs states. 
However, in the latter case, the value of tan /5 is assumed to be large, leading to strongly 
enhanced cross sections in the gg fusion and hh Higgs processes, while in the NMSSM one 
could have relatively moderate tan (3 values as in P5 and thus, smaller production rates. 

4.3 ATLAS strategy for NMSSM h\ a^al searches 

At the ATLAS collaboration [47], current efforts to find suitable search channels for special 
NMSSM phenomenological scenarios are being focused on the vector boson fusion produc- 
tion of a scalar Higgs boson with relatively low mass and subsequent decay via a pair of 
pseudoscalar Higgs bosons into four r-leptons, — > a^a^ — > 4r. This decay chain is favored 
in points P2 and P3 proposed here. The main emphasis is presently given to the case where 
all four r-leptons decay leptonically. 

A typical feature of the vector boson fusion production mode is the so-called tagging 
jets that are produced from the quarks that are scattered off the heavy vector bosons and 
merge to give the Higgs boson. These jets typically have high energies and lie in different 
hemispheres in the forward- and backward regions of the detector. Cutting on this signature 
is an important means to suppress background processes. Since there is no color flow between 
the quarks in the vector boson fusion process, jet production in the central detector region 
is suppressed. In contrast, central emission is favored in QCD interactions which constitute 
important background processes at the LHC [48]. Experimentally, this can be exploited by 
vetoing on additional jets in this region. 

The decay products of the Higgs boson typically lie in the central detector region. In 
general, leptons from the same pseudoscalar Higgs boson form pairs that lie close to each 
other in the detector, the separation being sensitive to the pseudoscalar mass. The invariant 
mass of the lepton pair has to be lower than the pseudoscalar mass, and is thus much lower 
than the Z^'-mass. For background processes including Z^-bosons, the photon interference 
therefore needs to be considered. In the experimentally most simple case, all four r-leptons 
decay to muons in the process a^a^ ^ 4r — ^ 4/i+4z/^+4z/T-. Since muons do not deposit 
considerable energy in the detector, muons from a very close pair can be also classifled as 
isolated. Decay channels including electrons need more consideration, as their energy is 
deposited in the calorimeters by electromagnetic showering. The possibility of separating 
nearby electrons or flnding un-isolated muons inside an electromagnetic shower from an 
electron needs careful study with a full ATLAS detector simulation. 

The transverse momentum of the stable leptons in this channel is rather low, since a large 
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part of the energy is carried away by the eight neutrinos in the final state. Therefore, not 
in all cases will all four leptons be identified. It might therefore prove favorable to require 
only three leptons to be found. For triggering, two muons (electrons) with px >10 GeV 
(20 GeV) or one muon (electron) with pt>20 GeV (25 GeV) are needed. It might also be 
considered to require a minimum transverse momentum for the remaining leptons to avoid 
having a high number of lepton fakes. Another feature of this channel is the large missing 
momentum from the eight neutrinos in the final state. 

In spite of the eight neutrinos in the final state and the fact that one lepton might 
remain unrecognized, it is still possible to reconstruct the mass of the scalar Higgs boson 
with help of the coUinear approximation [49] which is also used for mass reconstruction in the 
VBF, Higgs TT channel [50]. Since the pseudoscalar bosons and the r-leptons from their 
decays obtain a large Lorentz boost due to the large mass and pt of the Higgs boson, their 
decay products are emitted roughly in the direction of the original pseudoscalar. Exploiting 
momentum conservation in the transverse plane yields the 4-momentum vectors of the two 
pseudoscalars and thus the invariant mass of the scalar Higgs boson. The performance of 
this algorithm for the a^a? ^ 4r ^ 4/i + 8z/ channel is currently under study at 

ATLAS. 

Possible background processes for this channel are tt production, vector boson production 
in association with bottom or top quarks and production of vector boson pairs with additional 
light jets. Here, the leptons can come from decays of the heavy vector bosons or from decays 
of the bottom quarks. The tagging jets of the vector boson fusion channel might be faked by 
untagged 6-jets or by light jets from other sources in the event. It should be noted that the 
production of a vector boson pair in association with two light jets contains diagrams that 
have a structure similar to the vector boson fusion process, with two heavy bosons being 
scattered off by the incoming quarks and no color fiow between the quark lines. Possible 
methods to separate these background processes from the signal and their performance are 
currently studied by the ATLAS-collaboration. 

4.4 Prospects for the CMS experiment 

The final state topology — > a^a? t^t^t^t^ where the ~ 100 GeV h\ state is produced 
in the VBF process and the pseudoscalar a\ boson has a mass 2m,- < M^o < 2mh (so 
that the t^t~ decay mode is dominant) is currently under investigation by the CMS 

collaboration for the /^^/W^Tj^^rj^^ final state containing two same sign muons and two r jets. 
The r leptons from the decays of the light state are approximately coUinear and the 
non-isolated, di-muon high level trigger is needed to select the signal events. The standard 
CMS di-muon trigger with the relaxed isolation has the di-muon threshold of 10 GeV on 
both muons for a luminosity £ = 2 x 10^^ cm~^s~^ [52] and, since the muons from the signal 
events are very soft, the lower thresholds are needed. For instance, with a 7 GeV threshold 
the efficiency is increased approximately by a factor of two, but the QCD background rate 
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is also increased by the same factor [53] which is not acceptable. In order to cope with 
the rate, the same-sign relaxed di-muon trigger was recently introduced [54]. A PYTHIA 
simulation shows that the rate of di-muons from bb production is reduced by a factor of four 
when asking for the two same sign muons to have a threshold of 5 GeV. The off-line selection 
strategy requires the presence of the two same sign, non-isolated muons with one track in 
the cone around each muon direction, thus selecting the one prong r decays. 

For a point with masses M^o ~ 100 GeV and M^o ~ 5 GeV, one applies the following 
basic event selection cuts at the generation level: i) two same sign muons with px >7 GeV 
and |?7| <2.1 with one track of px >1 GeV in a cone 0.3 around each muon; ii) opposite 
charge for the muon and the track; in) two r jets with px >10 GeV and |?7| <2.1; iv) two jets 
with pt >30 GeV and |?7| <4.5. After these cuts, and using the SM VBF Higgs production 
cross section which is approximately 5.4 pb for the considered mass, the cross section of 
the signal after all selections is around 2 fb, leading to ~ 60 events at a luminosity of 30 
fb~^. The dominant backgrounds with two non isolated muons from the bb and tt production 
processes are under evaluation. 

For heavier bosons as in scenarios P2 and P3 where M^o ~ 10 GeV, the "non isola- 
tion" requirement should be relaxed since two r leptons from the a? — ^ tt decay are more 
separated and one should accept zero or one track in the cone (note that the relaxed di-muon 
trigger accepts both isolated and non isolated muons). One needs, however, to consider the 
backgrounds with isolated muons from W and Z° decays, in addition to bb and ti. 

The a\o!l — > t^t^t'^t~ decay with produced in Higgs-strahlung with leptonic 

decays of the W bosons, which can give a very clean and almost background free signal, is 
also being considered by the CMS collaboration. The leptons coming from W decays allow 
to trigger on the events and the CMS single isolated lepton trigger can be used [53]. The 
signal is unique as one has two coUinear r leptons due to the large boost of bosons and, 
potentially, the most interesting final state is when one of the rs decays hadronically and 
the other decays into muons. 

As for the selection criteria and the muon and r identification, for each r-jet candidate 
there must be a muon in the r-jet cone. The r jet and the muon are required to be oppositely 
charged and the events with two identified r jets with Et > 10 GeV and two muons with 
Pt > 7 GeV are selected. Unlike in the VBF case, the muons are not required to have 
the same sign. Finally, the events with extra jets on top of two r jets are rejected. A 
preliminary, full detector simulation and reconstruction analysis of the signal events in a 
scenario with M^^o ~ 100 GeV and M^o ~ 5 GeV shows that after all selection cuts, one 
expects ~ 10 events for 30 fb~^ data if the SM cross section of ~ 2.6 pb and a branching 
ratio BK{hi (^lO-i) — 0.9 are assumed. The potential backgrounds are ti,tW, QCD 
multijet, ly+jets and bb events. Preliminary tests with full detector simulation show that 
the requirement of two r jets, both having an oppositely charged muon in the jet cone, as 
well as jet veto suppresses the backgrounds very efficiently. 
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5. Conclusions 



The NMSSM is a very interesting supersymmetric extension of the SM, as it solves the 
notorious fi problem of the MSSM and it has less fine tuning. It also leads to an interesting 
collider phenomenology in some cases, in particular in the Higgs sector, which is extended 
to contain an additional CP-even and CP-odd state as compared to the MSSM. Hence, the 
searches for the NMSSM Higgs bosons will be rather challenging at the LHC in scenarios in 
which some neutral Higgs particles are very light, opening the possibility of dominant Higgs 
to Higgs decays, or when all Higgs bosons are relatively light but have reduced couplings to 
the electroweak gauge bosons and to the top quarks, compared to the SM Higgs case. These 
scenarios require much more detailed phenomenological studies and experimental simulations 
to make sure that at least one Higgs particle of the NMSSM will be observed at the LHC. 

In this note, we have proposed benchmark points in which these difficult scenarios are 
realized in a semi-unified NMSSM which involves a rather limited number of input param- 
eters at the grand unification scale and which fuUfils all present collider and cosmological 
constraints. In three of the five benchmark scenarios introduced here, the lightest CP-even 
Higgs boson decays mainly in two very light pseudoscalar Higgs states which subsequently 
decay into two b quarks or r leptons, leading to four fermion final states; in a fourth sce- 
nario, the next-to-lightest CP-even Higgs boson has a mass of ~ 120 GeV and couplings 
to fermions and gauge bosons that are SM-like but it decays into pairs of the lightest Higgs 
state, which then decays into a bb pair; in a last scenario, all neutral and charged Higgs 
particles are light (with masses less than ~ 180 GeV) and have weaker couplings to W/Z^ 
bosons than the SM Higgs particle. We have analysed the Higgs and supersymmetric particle 
spectra of these benchmark scenarios, discussed the various decay and production rates as 
well as other phenomenological implications and attempted to set the basis for the search 
strategies to be followed by the ATLAS and CMS collaborations in order to observe at least 
one of the neutral Higgs states in these scenarios. 

Addendum 

After the completion of this paper, discussions of various possibilities of the observation of 
the cascade decay aja? ... at hadron colliders appeared in Refs. [56]. 
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